Interference Study of the Xcoi^^Po) in the Reaction pp — ^ tt^tt" 

M. Andreotti,^ S. Bagnasco,^''^ W. Baldini,^ D. Bettoni,^ G. Borreani,'' A. Buzzo,^ R. Calabrese,^ 
R. Cester/ G. Cibmetto,^ P. Dalpiaz,^ G. Garzoglio,^ K. E. Gollwitzer,i M. Graham,^ M. Hu,i D. Joffe,^, 
J. Kasper,*^ G. Lasio/^'' M. Lo Vetere,^ E. Luppi,^ M. Macri,^ M. Mandelkern,'' F. Marchetto/ 
M. Marinelli,^ E. Menichetti/ Z. Metreveli,^ R. Mussa,^''^ M. Negrini,^ M. M. Obertino/'^ 
M. Pallavicim,3 N. PastroneJ C. Patrignani,^ T. K. Pedlar,^ S. Pordes,i E. Robutti,^ W. Roethel,6.4 
J. L. Rosen,^ P. Rumerio/'^ r. Rusack,^ A. Santroni,^ J. Schultz,* S. H. Seo,^ K. K. Seth,^ 
G. Stancari,i'2 m. Stancari,'*'^ ^ Tomaradzc,*^ I. Uman,^ T. Vidnovic III,^ S. Werkema^ and P. Zwcbcr^ 

(Fermilab E835 Collaboration) 

^Fermi National Accelerator Laboratory, Batavia, Illinois 60510 
^Istituto Nazionale di Fisica Nucleare and University of Ferrara, 44100 Ferrara, Italy 
^Istituto Nazionale di Fisica Nucleare and University of Genova, 16146 Genova, Italy 
"'University of California at Irvine, California 92697 
^University of Minnesota, Minneapolis, Minnesota 55455 
^Northwestern University, Evanston, Illinois, 60208 
''istituto Nazionale di Fisica Nucleare and University of Torino, 10125, Torino, Italy 

Fermilab experiment E835 has observed pp annihilation production of the charmonium state XcO 
and its subsequent decay into tt^tt". Although the resonant amplitude is an order of magnitude 
smaller than that of the non-resonant continuum production of 7r%o, an enhanced interference 
signal is evident. A partial wave expansion is used to extract physics parameters. The amplitudes 
J = and 2, of comparable strength, dominate the expansion. Both are accessed by L = 1 in the 
entrance pp channel. The product of the input and output branching fractions is determined to be 
B{pp Xco) X B{xcO vrV) = (5.09 ± 0.81 ± 0.25) x 10"^ 

PACS numbers: 13.25.Gv;13.75.Cs;14.40.Gx 



E835 studies charmonium formed in pp annihilation, 
a technique that allows direct access to all charmonium 
states. The experiment, located in the Antiproton Ac- 
cumulator at Fermilab, recorded new data in 2000. This 
run included 33 pb~^ of luminosity collected at 17 en- 
ergies (3340-3470 MeV) at the XcO- Results from the 
scan of the XcO resonance using the decay channel XcO ^ 
J/ip"f, J/ip — *■ e^e~ have already been published IJ. Si- 
multaneously recorded neutral events provide the tt'^tt'^ 
data reported here. The energy and luminosity measure- 
ments are the same for the J/V'7 and tt'^tt'^ channels. 

Charmonium resonances are scanned using a tunable, 
stochastically cooled p-beam that intersects a hydrogen 
gas jet target 9|. At each energy point, the cross section 
is measured by normalizing the number of events which 
satisfy the event selection criteria for the final state un- 
der study to the integrated luminosity collected. Typ- 
ically, an excitation curve is extracted from the large 
hadronic background by tagging electromagnetic final 
states. Hence, resonance parameters may be determined 
without relying on the detector resolution, but only on 
the knowledge of the p beam. At any setting of the 
beam momentum, the pp center of mass energy (Ecm) 
is known to 0.2 MeV; the spread (r.m.s.) in Ecm is on 
average about 0.4 MeV for these data. This resolution is 
substantially better than could be provided by detection 
equipment alone and is much smaller than the XcO width 



of 9.8 MeV m. However, when a large non-resonant con- 
tinuum is present in the final state of interest, as in the 
pp — > 7r°7r° analysis, special techniques must be used to 
determine the resonance parameters. 

In the vicinity of the XcO: the differential cross section 
for the process pp 



^OttO is 
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where Pj'^{z) are Legendre functions, x 
^^Xco)/rxco andz= Icosr 
production angle in the center of mass frame with respect 
to the p direction. 

The term ~Aii/{x + i) is the parameterization of a 
Breit-Wigner resonant amplitude with spin = 0. The 
partial wave sums in Eqs. (|2Jl and Q represent the con- 
tinuum contributions resulting from the initial states 
Ap — Apl = (helicity-0) and \Xp — Xp\ = 1 (helicity- 
1), respectively. These two helicity states are orthogonal 
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The Xco is only 



and do not interfere with one another, 
produced in the hehcity-0 initial state. 

At fixed z, the continuum terms Ae^^^ and B e^^'^ do 
not change markedly as the energy varies across the res- 
onance. It is then useful to rewrite Eq. 
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The above expression shows how even a small resonant 
contribution can lead to a detectable interference signal, 
albeit superimposed on a large continuum. For example, 
if the contribution Aj^ from the resonance at the reso- 
nance peak energy is 1% of the size of the cross section 
from the helicity-0 continuum A^, the contribution from 
the factor 2 An A will be 20 times larger than Aj^. The 
factor (sin (5^ — x cosSa)/ {x"^ + I) determines the shape 
of the interference pattern seen in the cross section. 

Two independent data analyses were performed and 
provide consistent results 0,3. One of them is presented 
here. The initial data sample consists of 4-photon events, 
selected by the neutral trigger (efficiency ^ 96.1%) and 
the central shower calorimeter. The energy resolution of 
this detector is as/E ~ 6%/./E(GeV) + 1.4%, while the 
polar and azimuthal angular resolution are erg ~ 6 mrad 
and 0-0 ~ 11 mrad, respectively A 5% cut on the 
confidence level of a 4C kinematic fit to the pp — > "f"f"f"f 
hypothesis ensures four-momentum conservation. Of the 
three possible ways to pair the 4 photons, one at most is 
a candidate for the pp — > tt'^tt'^ — > 7777 hypothesis due 
to the small opening angle of the 2-photon 7r° decays. 
A selection was made on the 2-photon invariant masses 
(to^j^^j and 771^3-),^) and the tt'^tt'' colinearity in the cm. 
frame. The total number of tt'^tt'^ candidate events is 
~ 500, 000. The product of the geometric acceptance 
and selection efficiency slowly increases up to z ~ 0.5 
(the average over this range is ~ 63%) where it starts 
to decrease rapidly. The instantaneous luminosity varied 
substantially during the data taking; the averages within 
each energy point were from 1.7 to 3 x 10^^ s~^cm~^. 
The event pileup was studied and corrected for by means 
of random triggers recorded throughout the data taking. 
The rate-dependent loss varied among the energy points 
from 14% to 23% with an average of ~ 20%. 

The background comes mostly from the 7r°7r°7r° and 
TT^cj — > 7r'^7r'^7 channels 0, and was determined by fitting 
the LEGO plot in Fig. [T| A --2% subtraction resulted 
from this background. 

The measured cross section of the pp tt^tt*^ reaction 
in the XcO region as a function of Ecm is shown in Fig.|21 
A non-resonant pp — > tt'^tt'^ production with a smooth 
energy dependence is present throughout the scanned re- 
gion. A clear resonance signal is visible close to the XcO 
mass; a finer scale plot shows a peak-shift of ~ 2 MeV 
(as a consequence of interference) toward low energy with 
respect to M^,,, = 3415.4 MeV/c^ of Ref. Q. 
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FIG. 1: The rrij-^j^ versus rrij^j^ LEGO plot in the region of 
the 7r°7r° peak, which is truncated at about 3% of its height. 
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FIG, 2; The pp tt^tt" cross section Act = /g^™"^ {da/dz)dz 
as a function of Ecm- The error bars are statistical. 



A binned maximum likelihood fit using the parame- 
terization of Eqs. 1^, lO and Q, with Jmax — 4, was 
performed simultaneously on all energy points. The XcO 
mass and width are set to the values (reproduced in Ta- 
ble that E835 measured via the reaction pp — > XcO 
J/7/'7, J/V' ^ e+e- Q, which had virtually zero back- 
ground and non-resonant cross section. The result of the 
fit is shown in Fig.O The curve A^+B'^ shows the sum of 
the non-resonant cross section contributions. The effect 
of the resonance, amplified by the interference, is seen in 
the separation (evident at small z) between da/dz and 
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FIG. 3: The pp tt' vr cross section versus z at Ecm ~ 
3415 MeV. The fit of Eqs. CJ, igj and ©, which was simul- 
taneously performed on all 17 energy settings, is shown along 
with its components. 
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FIG. 4: The pp tt^tt" cross section Act — J^^ {da / dz)dz 
as a function of Ecm- The reduced range fit and its compo- 
nents are also shown. 



+ and is almost entirely due to the interference- 
term of Eq. 101). The separation decreases as z increases, 
following the trend of A. The term is small at small 
values of z, due to a factor z present in all the associate 
functions Pj(z). The net suppression factor of B^ with 
respect to A^ is z^ at small z. The contribution of the 
"pure" resonance, Aj^/{x'^ + l)i is negligible and not dis- 
tinguishable from zero in the figure. 

The fit has 408 bins: 17 energy points times 24 bins in 
z (from to 0.6). The number of free parameters is 15: 
the resonance amplitude, Afi; the coefficients C,/=o.2,4 
and Cj^2 4 (each of them is given a linear energy depen- 
dence); and the phases (5j=o,2,4 and {Si — (^2)- Including 
energy dependence in the phases does not significantly 
alter nor improve the fit. 

The partial wave summation is truncated at Jmax = 4. 
The addition of a relatively small J = 4 amplitude 
to a basic Jmax = 2 fit is already "fine tuning" . At 
Ecm = 3415 MeV the relative amounts are Cq : (C2, C2 ) : 
(C4,Ci) ~ 1 : (0.65,0.41) : (0.20,0.12). J = 6 and 
higher partial waves are highly oscillatory functions and 
do not significantly improve the fit nor are physics moti- 
vated. Lpp — 1 "feeds" J = and 2. Lpp = 3 is the min- 
imum entrant angular momentum for J = 4. Increasing 
Lpp correlates directly with the collision impact parame- 
ter (5). The physics of the annihilation process into tt^tt^ 
at small z strongly favors small b values. The Xco channel 
obviously requires total valence quark annihilation. The 
non-resonant tt^tt*^ channel does not require total valence 



quark annihilation; however, hard, short-range collisions 
are required to redirect the non-annihilating quarks into 
reforming as part of the tt'^tt'^ final state. Our truncating 
partial wave expansion supports this physical picture. 

In order to extract the product of the branching ra- 
tios B{xcO pp) X B{XcO — > 7r°7r°), the natural place to 
focus on is the small z region, where the suppression 
of the non- interfering helicity-1 continuum (of magnitude 
B^, see Eq. ) with respect to the interfering helicity-0 
continuum (of magnitude A^) is exploited. At small z 
the uncertainty on the relative amounts of A^ and B^, 
critical in quantifying the amplification effect of the in- 
terference, is severely limited and guarantees a model- 
insensitive result. In addition, the interference-enhanced 
XcO signal has a substantial size at small z and decreases 
at increasing 2, while the forward-peaked continuum pro- 
duction dominates at larger z. Accordingly, a new fit was 
performed in a reduced range < z < 0.125. Eq. 
was used again with the following differences. First, the 
small contribution of B^ — \Be^^^ p was fixed and taken 
from the fit in Fig.|31 which is dominated by the forward 
peak of this component. Second, a new parameterization 
of the helicity-0 continuum was introduced; due to the 
small range of z, the number of free parameters can be 
decreased by employing a polynomial expansion on z and 
X for = ao + aix + a2X^ -|- a^z^ . This is an adequate 
approximation to the partial wave expansion in the small 
z interval. The resulting da/dz, integrated over the re- 
duced range, is shown in Fig. ^ As evident from the 
figure, the constant term oq is dominant, while ai and 
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TABLE I: E835 Results (errors are statistical and systematic, 
respectively) . 



Bin — 

Bout = 


Common cham 

Bixco-^ J/i-jT 


lel B{xco -» pp) 
B{x.o ^ ^V«) 


M^^, (MeV/c^) 
Txco (MeV) 
X Bout (10-^) 


3415.4 ±0.4 ±0.2" 
9.8 ±1.0 ±0.1" 
27.2 ± 1.9 ± 1.3* 


3414.7+og ±0.2" 
8.611J ±o.r 

5.4:2t°oil ± 0.22"= 


Final result for Bi„ x Bout (10 ') 
and phase Sa (degree) 


5.09 ±0.81 ±0.25" 

39 ± 5 ± 6"^ 



"The J hp was detected through its decay into e+e QQ. 
"From where Bi„ xBout, ^Xco l^Xco were free parameters. 
'^This analysis with Bin xBout, J^xca ^Xco ^ ^""^^ parameters. 
''This analysis with M^^q and r;^^Q set to values from Ref. 



a2 provide a small slope and curvature as a function of 
X. The small dependence in the reduced z region is 
accommodated by 03. The fit has 6 free parameters (A^?, 
ao, ai, 02, as, and Sa) for 85 bins (17 energy points times 
5 bins in z, from to 0.125). By searching for improve- 
ments in the it is found that the phase 5a does not 
exhibit any dependence on z in this small 2-range nor 
on the energy, and that does not require additional 
powers of x and z, nor mixed terms such as xz'^ and x'^z^. 

It has been noted that the "pure" Breit-Wigner (the 
fictional cross section that would result if the non- 
resonant amplitudes could be turned off) is very small. In 
Fig. 21 it is shown enhanced by 20 to make it comparable 
to the signal actually detected. 

The fit provides the value of An, which is related to 
the product of the input and output branching ratios by 
A% = nX^ X B{xcO ^ pp) X B{xcO Tr^Tr"), where A 
is the center of mass de Broglie wavelength of the initial 
state. The result for the product of the input and out- 
put branching ratios obtained with M^^^ and T^^q con- 
strained to the values obtained in our J/iJj 7 measurement 
Q is reported in TableHjas "Final result for Bin x Bout" ■ 
The dominant systematic error arises from the luminos- 
ity determination. The uncertainty on the knowledge of 
the helicity-1 continuum affects Bin x i?o«t by ~ 1%. 

The difference of phase (Sa) between the amplitudes 
of the helicity-0 continuum and the resonance produces a 
modest constructive interference on the low-energy and 
destructive on the high-energy side of the Xco mass. 

A consistency check is provided by allowing M^^^ and 
^Xco to vary along with Bin x Bout and 6a- TableQ]shows 
that the resultant M^^^g and F^^^ values agree with the 
J/ipj values. Although the tt^tt*^ sample is more copious 
than the J ftp"/, the values determined by the ttOttO fit 
have larger uncertainties because of the higher number 
of (coupled) fit parameters. 

The results presented so far are obtained from the E835 
2000 data sample alone, with the exception of the well- 
known B{J/tp — > e+e~) 0. To deconstruct the entrance 
and exit channel branching fractions, we must use data 



from the literature. B{xcQ — > tt^tt") is more than 
an order of magnitude larger (and thus better measured) 
than B{xcO — > Pp)- We then obtain B{xcO —>■ Pp) = 
(2.04 ± 0.32,t,t ± QAQsyst ± 0-2Spdg) x IQ-*. 

In addition, taking the ratio of Bin x Bout f^or 
the two channels we measured, we obtain B{xcO — > 
J/^-f)/B{xcO TfOTfO) = 5.34 ± 0.93 ± 0.34 (a num- 
ber of minor systematic uncertainties cancel). Using 
we then determine B{xcO Jli^l) — (1.34 ± 
0.23stat ± 0.09syst ± 0.19p_dg)% and, taking T^^„ from 
Q, r^^„_j/^^ = (131 ± 26«tat ± isyst ± 18pdg) kcV. It 
is interesting to compare the El radiative transition of all 
three Xcj states. The above measurement of T^^^^^j/^^ 
is now in excellent agreement with the energy indepen- 
dent scaled quantities T^^j^j/^^/qj for Xci and Xc2 0- 

Summarizing, an interference pattern in the tf^tf" cross 
section has been observed at the XcO mass. An origi- 
nal analysis has been developed to extract the XcO res- 
onance parameters. Combining the present study with 
our previous one of — > J/V' 7 Q , important improve- 
ments in the knowledge of the XcO are achieved. The 
presented work proves that resonances can be observed 
and studied via interference in final states dominated by 
non-resonant channels. The developed analysis could be 
adopted in future studies, such as tf^tf" and tt^t] scans to 
search for possible DD bound systems at ~ 3700 MeV. 
Other two-body final states could be employed for the 
study of charmonium singlet states via interference. 
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